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Abstract 

Configuration-Interaction (CI) models are approximations to the 
electronic Schrodinger equation which are widely used for numerical 
(— I ' electronic structure calculations in quantum chemistry. Based on our 

0^- recent closed-form asymptotic results for the full atomic Schrodinger 

r~| ! equation in the limit of fixed electron number and large nuclear charge 

|FG09j , we introduce a class of CI models for atoms which reproduce, 
at fixed finite model dimension, the correct Schrodinger eigenvalues 
and eigenstates in this limit. 

We solve exactly the ensuing minimal model for the second period 
^ ' atoms, Li to Ne. The energy levels and eigenstates are in remarkably 

, good agreement with experimental data (comparable to that of much 

larger scale numerical simulations in the literature), and facilitate a 
mathematical understanding of various spectral, chemical and physical 
If^ . properties of small atoms. 

o 

^ : 1 Introduction 

> . 

^ ' From the early days of quantum mechanics it has been clear that the chemical 

. behaviour of atoms and molecules is governed by their energy levels and 

electron configurations, which in turn are determined, to very high accuracy, 
by the eigenvalues and eigenstates of the Schrodinger equation Hip = Eip. 
But 80 years on, high-accuracy numerical computation of such data remains 
a largely unresolved challenge, even for the smallest of systems such as a 
single Carbon atom. The only computations of which we are aware which 
meet the mathematical ideal [BLWW04] of convergence tables showing an 
increasing number of converged digits as a function of basis set size or number 
of iteration steps (for a reproducibly documented algorithm for the original 
problem) concern two-electron systems such as He and H2. See |KNN08j for 
recent advances and references. 

The underlying reasons are two-fold. 
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First, a "curse of dimension" phenomenon is present: the Schrodinger 
equation for an atom or molecule with electrons is a partial differential 
equation in R^^, so direct discretization of each coordinate direction into K 
gridpoints yields K^^ gridpoints. Thus the Schrodinger equation for a single 
Carbon atom (A^ = 6) on a ten point grid in each direction {K = 10) already 
has a prohibitive 10^^ degrees of freedom. 

Second, one is dealing with a tough multiscale problem: chemical be- 
haviour is not governed by total energies, but by small energy differences 
between competing states. Even for very small systems, these are typically 
several orders of magnitude smaller than total energies. For instance, as 
shown in the table below, the spectral gap between ground state and first 
excited state of the second period atoms is less than 1 % of the total size of 
these energy levels in all cases, and only about 0.1 % for Carbon, Nitrogen 
and Oxygen. Nevertheless this tiny gap is of crucial importance, as the two 
states it separates have different spin and angular momentum symmetry, and 
hence completely different chemical behaviour. 



Atom 


Li 


Be 


B 


C 


N 





F 


Nc 


Ratio of first spectral gap 
to ground state energy 


0.0093 


0.0068 


0.0053 


0.0012 


0.0016 


0.00096 


0.0078 


0.0047 



Table 1: Multiscale structure of atomic spectra (experiment, |RJK"'"07] ). 



To deal with the curse of dimension, in quantum chemistry a large array 
of reduced models has been developed. For small systems with up to one or 
two dozen electrons, the most accurate and most widely used class of models 
are the Configuration-Interaction (CI) models, whose origins go back to the 
early years of quantum mechanics (see e.g. |Hyl29| ) and whose systematic 
development started with the work of Boys |Boy50| and Lowdin |Loe55j . 
Roughly speaking, these are "tensor product Galerkin approximations" : the 
full electronic Schrodinger equation is projected onto a subspace spanned by 
carefully chosen Slater determinants (= antisymmetrized tensor products), 
which are in turn formed from a small set of orbitals (= elements of the 
single-particle Hilbert space L^(]R^ x Z2)). 

Different CI models differ by the choice of orbitals and the selection of 
the subset of Slater determinants. The question of how to best make these 
choices remains the subject of a great deal of current research in the quantum 
chemistry literature, with the best methods to date relying on a combination 
of chemical intuition, computational experience, and nonlinear parameter 
optimization, as well as on a huge number (between 10^ and 10^) of included 
determinants. See [SU96[ IHJOOOj for a general overview of the CI method 
and its most common variants such as Doubly excited CI (DCI), Multi- 
determinant Hartree-Fock (MDHF), Complete active space self-consistent 
field method (CASSCF), Coupled- Cluster theory (CC), and the (desirable 
but usually not practical) Full CI (FCI), and see e.g. |BT86[ nT5T94l IKR021 
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IBMn4t K^NMCm^t Uoh.05\ INNKIOTl IKNNOS] for applications to atomic en- 



ergy level calculations. 

Our goal in this paper is to introduce, analyze, and apply to atomic energy 
level prediction a particular class of CI models for atoms which exploit our 
recent closed-form asymptotic results for the full atomic Schrodinger equation 
in the limit of fixed electron number and large nuclear charge Z |FG09] . 
Namely, we require that the model of fixed finite subspace dimension K 
reproduce correctly the first K Schrodinger eigenvalues and eigenstates in 
this limit. 

That such a requirement can be met by a fixed-resolution CI model is not 
trivial (for example, it is not met by Hartree-Fock theory, even in an infinite, 
complete one-electron basis), but a simple consequence of the asymptotic 
results in |F(;n9j (see Section [23]) . 

The above limit exhibits the important multiscale effect that the ratio 
shown experimentally in Table [T] of first spectral gap AE to ground state 
energy E of the Schrodinger equation tends to zero |FG09j . The require- 
ment that the corresponding eigenstates and gaps be nevertheless captured 
correctly by an approximation should hence be relevant to yielding good 
eigenstates and gaps in the realistic situation when this ratio is small. 

In fact, even the minimal asymptotically correct CI model for atoms and 
ions with 1 to 10 electrons (eq. (A'), (B'), (C) in Section (2111), whose subspace 
dimension turns out to be 8, 28, 56, 70, 56, 28, 8 for Li, Be, B, C, N, F, turns 
out to be very interesting. 

(i) We find that the requirement of asymptotic correctness leads to Slater 
orbitals P(x)e~'^'^' (where P is a polynomial and 7 > a constant), 
not Gaussian orbitals P(x)e~"'^^^ used in the overwhelming majority 
of numerical CI computations on account of their easy facilitation of 
two-centre integral evaluation. See Section 12.41 

(ii) The model is exactly soluble (via the methods introduced in |FG09] ). 
except for the determination of the dilation parameters 7 in the Slater 
orbitals, which are straightforward to calculate numerically with very 
high accuracy. See Section [3l 

(iii) The model does remarkably well when compared to experimental data 
and high-dimensional simulations in the literature. It captures around 
99 percent of the ground state energy in all cases, without a single 
empirical parameter! See Figure [H Moreover the predicted ground 
state spin and angular momentum quantum numbers (^S for He, Be, 
Ne, 2S for H and Li, ^S for N, for B and F, and for C and O) 
come out right in each case; spectral gaps are captured well (and in 
a significant number of cases more accurately than in the benchmark 
numerical multi-determinant Hartree-Fock calculations of Tatewaki et 
al |TTST94] which used a much larger basis set); and for > 5 the 
model is never outperformed by more than a factor ten by any method. 
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including large-scale simulations with subspace dimension bigger than 
10''. For a detailed comparison see Section HI 



Atom 

He Li Be B C NO F Ne 




Figure 1: Experimental (dashed) and minimal asymptotics-based CI (solid) 
ground state energies of atoms. 

Thus, our work yields for the first time few-parameter, explicit, closed-form 
approximations to the low-lying eigenstates of the atoms Li, Be, B, C, N, O, 
F, Ne which are of chemically relevant accuracy. These provide a hopefully 
useful reference for the calibration of numerical methods, and a valuable tool 
to advance mathematical understanding of physical, chemical and spectral 
differences between the elements. For example, the ground state wavefunc- 
tions confirm the basic mathematical picture of the periodic table obtained 
in |FG09] by asymptotic analysis of the Schrodinger equation for strongly 
positive ions, and make it quantitative for neutral atoms; and they allow to 
trace the size of spectral gaps to individual Coulomb and exchange integrals, 
thereby making more rigorous longstanding insights by quantum chemists 
and revealing the cancellations that lead to the small size of gaps compared 
to total energies (Table [1]). See Section HI 

Nevertheless a great many open problems remain, even for minimal asymp- 
totics-based CI. 

1) In this paper we demonstrate its accuracy via comparing to exper- 
imental values (see Section Hj) and proving desirable theoretical properties 
(see Section [2]), but how can it be understood in terms of rigorous error 
estimates comparing it to the Schrodinger equation? 

2) In particular, why is the use of just one dilation parameter per one- 
electron subspace Vj so effective? As far as we are aware, although screening 
parameters are widely used (see Section [2]), there are no rigorous mathemat- 
ical results regarding their effectiveness. For instance, one might hope such 
parameters to emerge in some order expansion of a suitably scaled problem. 

3) How does the model fare for larger atoms? For this step we would sug- 
gest automation of the calculation of the eigenspaces and energy expressions 
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(analogous to Tables [2H3]) , and Fourier transforms and one- and two-body 
integrals (as in Lemmas 13.11 and I3.2p . We would hope the model to show 
interesting chemical effects such as the shell ordering 4s < 3d, and its occa- 
sional reversal, in the transition metals. 

Finally, it is highly desirable that the asymptotics-based CI approach in- 
troduced here be extended to molecular problems. The principal observation 
(Theorem 12.11 (ii)) that CI models of fixed finite subspace dimension can be 
constructed which reproduce correctly the first K Schrodinger eigenvalues 
in a large nuclear charge limit is not limited to atoms, as will be discussed 
elsewhere. But in the molecular case the ensuing orbitals are not available 
in closed form, and hence do not lead so readily to a mathematical picture 
of basic physical and chemical properties. 

2 Asymptotics-based CI models for atoms 
2.1 General CI models 

We begin with a mathematical description of CI methods. We find it con- 
venient to do so in the more abstract setting of subspaces and subspace 
projections rather than the, equivalent, setting of basis sets and expansion 
coefficients used in the chemistry literature |S096j . Moreover we introduce 
a rigorous distinction between general and symmetry-preserving CI meth- 
ods. (Both of these, as well as hybrid methods in which the solution to a 
non-symmetry-preserving model is projected a posteriori onto an invariant 
subspace, are in use in the chemistry literature.) 

Starting point for the derivation of any CI model is the exact (non- 
relativistic, Born-Oppenheimer) time-independent Schrodinger equation for 
atoms and ions which one seeks to approximate, 

i/^ = E^, (1) 

where, for nuclear charge Z > and electrons and in atomic units, 

i=l ' l<i<7<Af ' * 
^ ^ ' V ' 

E eR, and 

<f ELldR^xZ^f). (3) 

Here and below the Xi G M.^ are electronic position coordinates, Sj G Z2 = 
{±|} are spin coordinates, and is the usual Hilbert space of A^-electron 
functions \E' : (M^ x 1^2)^ C which are square-integrable, 

/ V |^(xi,si,...,a;iv,sjv)r = ll^'ll^ < 00, (4) 
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and satisfy the antisymmetry principle that, for all i and j, 



. . y Xi, . . . , Xj , s J , . . . ) ( . . . , Xj , Sj , . . . , Xi, Si, . . . ^ . (5) 

Mathematically, if is a bounded below, self-adjoint operator with domain 
L'^nH'^, where is the usual Sobolev space of functions with second weak 
derivatives belonging to L^. It is known that for neutral atoms {Z = N) and 
positive ions {Z > N), there exists an infinite number of discrete eigenvalues, 
the corresponding eigenspaces being finite-dimensional (Zhislin's theorem, 
see |Fri03j for a short proof). 

Translating }S096j into mathematical terminology, a CI model is a tensor 
product Galerkin approximation to the many-electron Schrodinger equation. 
More precisely: 

Definition 2.1 A CI model of an A^-electron system with Hamiltonian H is 
a projection of the Schrodinger equation ([T]) of form 

PHP^i = E*, ^ G V, V a subspace of Ll, P = projector onto V, (6) 

with the additional requirement that V must possess a basis consisting of 
Slater determinants. 

Recall that a Slater determinant is an anti-symmetrized tensor product \rji . . . t/n) G 
L^((R^ X 7^2)^) of orthonormal one-electron functions rji G L^(M^ x Z2), the 
antisymmetrization being necessary to comply with the quantum mechanical 
law i^. The difference between different CI models lies in the freedom to 
choose the subspace V, or - in quantum chemistry language - to select a 
set of orbitals and a set of Slater determinants to be included into the CI 
expansion. 

Note that if V is spanned by the orthonormal Slater determinants \E'a, a = 
1, . . . , M, the projection operator P onto V has the expansion "^^^i \'ifa){'^a\, 
and eq. (E]) can be written in its more standard matrix form Ac = Ec, where 
A is the MxM matrix with entries Aajj = {^alHl"^ f^) , and c is the coeffi- 
cient vector in the expansion = Ca'^a- The more abstract form ([6]) 
emphasizes the elementary fact that the CI eigenvalues and eigenstates only 
depend on the subspace V, not on the choice of basis within this subspace. 

A basic desirable feature of CI models, not related to the tensor product 
structure but only to that of a linear subspace projection, is the following. 

Lemma 2.1. (Rigorous upper bound) For j = 1, . . . , dimV, we have Ej < 
Ej'^ , where Ej and E^^ are the eigenvalue, including multiplicity, of the 
Schrodinger equation ^ respectively the CI equation 

Proof This is an immediate consequence of the min-max theorem for dis- 
crete eigenvalues of a self-adjoint operator below the bottom of the essential 
spectrum |RS78j . 
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2.2 A mathematical definition of the notion of config- 
uration for atoms 



In the quantum chemistry hterature, the word "configuration" is often em- 
ployed as a synonym for Slater determinant |SQ96j . But in the atomic 



physics and atomic spectroscopy literature (e.g. [RJK+O?] ). as well as some 
of the best computational studies, the word "configuration" has a more sub- 
tle meaning, which takes into account the important role played by spin and 
angular momentum symmetries. For our mathematical purposes, the latter 
notion turns out to be very useful, so let us formalize it mathematically. 

First, recall the total angular momentum operator L = (Li, 1/2,-^3), the 
total spin operator S = {Si, S2, S^) and the parity operator R, along with 
the fact that the operators 

L^S^L3,^3 and ^ (7) 

commute with each other and with H (see |FG09j for the result, as well as a 
mathematical definition of the operators L, S and (171)) Q 

One starts from a finite number of mutually orthogonal subspaces of the 
single-electron Hilbert space, 

Vi,...,VkCL\R^xZ2), (8) 

which are irreducible representation spaces for the joint spin and angular 
momentum algebra Span {Li, L2, L^, Si, S2, S^}. In elementary terms, this 
means that the subspaces must be of "fixed angular and spin symmetry" 
and "minimal dimension", more precisely: each Vi must be invariant under 
the Li and Si, the operators and must be constant on Vi, and Vi must 
have minimal dimension (i.e. dimension (2s + + l)|s=i/2 = 2 ■ {2i + 1) 
when = 1)1 and = s{s + 1)1; note that the spin quantum number 
s equals 1/2 for any Vi, since = |/ on the whole single-electron state space 
L2(R3 X Z2)). 

Definition 2.2 A configuration of an A^-electron atom or ion is a subspace 
of N-electron state space of the following form: 

Cdi,-,d, ^ Span||xi,...,X7v) {xi,---,Xn} any ON set with ^{i\xi G Vj} = rf^}, 

(9) 

where Vi, .., C L^(]R'^ x Z2) are mutually orthogonal irreducible represen- 
tation spaces of the joint spin and angular momentum algebra, and di, ...,dk 
is a partition of (i.e. rfj G N U {0}, Y.j dj = N). 



^On single-electron functions xi^^s), x G M"^, s £ {±1/2}, one has L = a; A iV, 
{Rx){x, s) = x(— a;, s), and Sa is multiplication by a Pauli matrix, 

/ {So.x){x,l/2) ( x{x,l/2) \ ^ 1 /O 1\ ^ 1 /O ~i\ ^ 1 A 



(^„x)(^,-i/2);"^nx(^,-i/2);''''^~2U 0i'^2^-2U 2V0 -1, 
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The main point here is that all choices of the x^'s consistent with the re- 
quirement that a fixed number of them have to be picked from each Vj have 
to be included. 

As an elementary but important consequence, each configuration is in- 
variant under the spin and angular momentum operators L and S, and 
in particular under the operators (171). This is immediate from the invari- 
ance of the Vi under L and S and the following identity for the applica- 
tion of one-body operators to Slater determinants: ' ' ' Xn) = 

J2j=i\Xi---{Axj)---XN)- 

Example 1: The configurations ls^2s and ls^2p for Lithium. 

Let 

Vi = Span {01s T, 0i4}, V2 = Span {024, 024}, 

Vs = Span {02p, T, 02p, i M = 1, 2, 3}, (10) 

where the 0's are the hydrogen-like orbitals 

73/2 

0is(x) := ^e-^^N 
02s(a;) 



(mIz?^^)''' (1 - 1(2^1 + Z2)\x\) e-^^N/^ (11) 
^5/2 

02p,(a;) := -^x,e~^-^\^\/', i = 1,2,3, 

Zi, Z2, Zj, are positive parameters, and t, | denote the spin functions { (s) = 
5i(.), i(.)=5_.(.). 

Note that the orbitals in (fTOj) are orthonormal (hence the coefficients in 
02s), and that for Zi = Z2 = ^3 = Z they reduce to the standard eigenstates 
of the hydrogen atom Hamiltonian — |A — y|| . The Zi will play an important 
role later. 

For = 3, choosing the partitions di = 2, ^2 = 1, c^s = respectively 
(ii = 2, 6/2 = 0, (is = 1 yields the subspaces (or configurations) 

C^'l'O = Span{|0i4, 0i4, 02s T), I0l4, 0l4, 02s i)}, 

C2'0'i= Span{|0i4, 0i4, 02pJ), 101s T, 0isi,02p4) M = 1,2,3}. (12) 

We call these subspaces ls^2s^ and ls^2p^. In chemistry this terminology 
is common to describe the structure of individual wavefunctions, but in the 
setting just introduced, it is independent of which wavefunction is chosen. 

For ^1 = ^2 = ^3 = Z, these subspaces have the interesting physical 
meaning that they are the bottom two eigenspaces of the Lithium atom 
Hamiltonian in first order perturbation theory [FG09j . 

Example 2: The subspace 2p\2pl for Helium 
The subspace 

C= Span{|02p4,02p4), l02p4,02p2T), l02p4,02p2T), I02p4,02p2i)} 
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is not cL configuration, because the selection of Slater determinants does not 
correspond to the rule in Definition 2.1. Indeed this subspace is not invariant 
under the spin and angular momentum algebra. For instance, applying L3 
to the first Slater determinant gives i(|02pi T 02pi i) — |02p2 T 02p2 i))) which 
lies outside the subspace. 

2.3 Symmetry-preserving CI models 

A general class of symmetry-preserving CI models can now be defined math- 
ematically. We remark that the principle of symmetry-preserving numerical 
schemes has proved very successful in other areas of scientific computing, a 
prime example being symplectic schemes in Hamiltonian dynamics |LR05] . 

Definition 2.3 A symmetry-preserving CI model for an A^-electron atom or 
ion with Hamiltonian H is a finite-dimensional projection of the Schrodinger 
equation ([1]), 

PHP^ = E"^, \l/ G V, V a subspace of L^, P=orthogonal projector onto V, 

(13) 

with the additional requirement that 

= Span{C«,...,C(^)}, (14) 

where Vi,..,Vk C L^(]R^ x Z2) is a collection of mutually orthogonal irre- 
ducible representation spaces of the spin and angular momentum algebra, 
and each C*^-'^ is a configuration with respect to the Vi. 

Example Taking V = Span {C^'^'° , C^'O'i} = Span{ls22s, ls^2p} (with no- 
tation as in Example 1) yields an invariant CI model for Lithium. 

The fundamental point of Definition 2.3 is that unlike general CI, symmetry- 
preserving CI retains the spin and angular momentum symmetries of the 
atomic Schrodinger equation. In particular, eigenspaces retain well defined 
spin and angular momentum quantum numbers 5" and L (see |FG09j for their 
mathematical definition) : 

Lemma 2.2. (Symmetry preservation) For arbitrary N and Z, and any 
invariant CI model [T^) . [14\ ) as in Definition 2. 3, The operators ^ 
(i) leave V invariant 

(a) commute with the CI Hamiltonian PHP : V ^ V . 

In particular, any eigenspace of an invariant CI model possesses a basis of 

joint eigenstates of the operators 

Proof This is an elementary consequence of the invariance of individual 
configurations under L, S and R. Note that the underlying one-electron 
subspaces Vi, being L^ eigenspaces with some eigenvalue i{i + 1), are auto- 
matically R eigenspaces, with eigenvalue (—1)^"^. 

We note the well known fact that the physically important property (i) is vi- 
olated by standard approximations such as the Hartree-Fock approximation. 
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even when the individual orbitals have well defined spin and angular momen- 
tum quantum numbers. For instance, the Slater determinant |02pi T 4>2p2 i) 
is neither an nor an eigenstate. 

2.4 Asymptotics-based subspace selection 

We now come to the, in applications crucial, issue of selecting a "good" CI 
subspace V in the approximation (Q. 

Commonly, this relies on a great amount of chemical intuition, computa- 
tional experience, and nonlinear optimization. For example, one would em- 
ploy the set of Slater determinants formed from the first N + k eigenstates of 
the nonlinear Hartree-Fock equations of the system under consideration ( "k- 
fold excited CI" ) , solved numerically in a background subspace of dimension 
L » N + k spanned by Gaussian orbitals. For more information, common 
variants and refinements see |S096t IHJOOO] . 

We propose here an alternative strategy, in which the intermediate step 
of a Hartree-Fock calculation no longer appears, and which is based on three 
reasonable theoretical requirements. The CI model should 

1. preserve the symmetry of the atomic Schodinger equation under spatial 
and spin rotation (see Definition 2.3) 

2. preserve the virial theorem, i.e. eigenstates should have the correct 
virial ratio of —2 between potential and kinetic energy 

3. be asymptotically correct in the iso-electronic limit Z ^ oo. 

By 3. we mean that the model (if its dimension is K) reproduces correctly the 
first K Schrodinger eigenvalues and eigenstates in this limit (see Theorem 12. II 
for a precise statement). Note that the limit of large Z captures the physical 
environment of inner shell electrons in large atoms. Also, recall its important 
theoretical feature that the ratio of first spectral gap to ground state energy 
of the Schrodinger equation tends to zero |FG09j . the experimental 

ratio for true atoms being very close to zero (see Table [1]). 

We now apply requirements 1, 2, 3 to the atoms Li to Ne, by not designing 
a largest such model which can be handled computationally, but a minimum- 
dimensional model. Below, {Zi, Z2, Z^) =: Z_ denotes the vector of dilation 
parameters appearing in the orbitals (|TT|) . We first discuss the case of the 
ground state. The ensuing minimal CI model for Li, Be, B, C, N, O, 
F, Ne ground states is then: 

(In (A), it is understood that only configurations C'^^''^2,di ^qj- which each 
cij is > are included.) 

Some remarks are in order. 

(1) This model is certainly not the only conceivable model which satisfies 
1., 2., 3., especially since condition 3. is only asymptotic, but it is probably 
the simplest. The subspace V{Z^) in (A) comes from the theorem in [FG09j 
that the above subspace with Zi = Z2 = Z^ = Z is asymptotically equal to 
the union of the lowest eigenspaces of the full Schrodinger equation ([T]). In 
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(A) (Choice of a parametrized, asymptotically exact family of subspaces) 
V{Z) := Span{C2.0'^-2,C2'i'^-3,C2'2.^-4}, 

with the C's as in (El) and Vi, F2, V^s as in ((lOD, ([U]), 

(B) (Subspace eigenvalue problem) 

E^\Z) := lowest eigenvalue of P{Z)HP{Z) on F(Z), 

\I/'-^-^(Z) := corresponding normalized eigenstate, 

where P{Z) = orthogonal projector of L^^^iW" x 1.2)^) onto F(Z) 

(C) (Variational parameter determination) 

:= argmin2E^^(Z), E^' := E^'iZ,), ,= ^ci(^zj 



particular, this theorem dictates that the Vi should consist of Slater orbitals, 
not the commonly used Gaussian orbitals. The presence of the variable 
dilation parameters Zj and eq. (C) comes from requirement 2., which is 
equivalent to stationarity of the energy of eigenstates with respect 

to dilations (see the proof of Theorem 12. ip . 

(2) There is no empirical parameter. 

(3) The model has the following variational formulation: 

E^^ = min min 

ze(o,oo)3 *ey(z), 11*11=1 

with the set of minimizers \l/ being equal to the set of normalized lowest 
eigenstates of (B). This is an immediate consequence of (C) and the Rayleigh- 
Ritz variational principle for the bottom eigenvalue in (B). 

(4) Dilation parameters like the Zi are closely related to physical ideas 
of screening, and go back at least to Slater (in the context of the Hartree 
equations, [SlaSOj ISla64j ). They are widely used in the quantum chemistry 
literature, and are in most studies determined a priori, e.g. via a Hartree- Fock 
calculation (see |BT86t IS096j ). However, from a mathematical standpoint 
it is of interest to determine them variationally for each eigenstate, as done 
here; this implies that the ensuing wavefunctions satisfy the virial theorem 
(see Theorem 12. ip . Note also that validity of the latter cannot be guaran- 
teed by linear parameters (i.e., subspace enlargement), but requires making 
the model nonlinear. This is because the dilation group 

which underlies the virial theorem, is a non-compact group which - unlike 
the compact groups 5*0(3) and SU{2) corresponding to angular momentum 
and spin - leaves no finite-dimensional subspace of L^^j-((M^ x Z2)) invariant. 
The proof of this fact is left to the interested reader. 

We now extend (B), (C) to excited states. The simplest generalization would 
be to compute all eigenvalues Ej{Z_) and corresponding orthonormal eigen- 
states ^I/j(Z) of P{Z_)HP{Z_), then minimize each eigenvalue over Z_. 

But this procedure does not maintain the basic property of the full Hamil- 
tonian that eigenstates with different eigenvalue are orthogonal. However, 
the symmetries described in Lemma [2.21 come to our help. If two eigenstates 
of the CI Hamiltonian are also simultaneous eigenstates of the operators ([7]), 
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which we can assume by Lemma 12.21 then they remain orthogonal after min- 
imization of their eigenvalues over the Zj, as long as the eigenvalue of at least 
one of the operators ([7]) are different. Thus, in each symmetry subspace (i.e., 
each joint eigenstate of the operators ([7])) we determine the values of the Zi 
that yield the minimum value for the lowest eigenvalue in the subspace, then 
use this value to calculate all eigenvalues and eigenstates in the subspace. 
This way, orthogonality is maintained and in particular the CI energy lev- 
els remain rigorous upper bounds to the true energy levels. In practice this 
method is very close to minimization of each eigenvalue, since most sym- 
metry subspaces turn out to be one-dimensional, and none are more than 
two-dimensional (see the next section). The use of the Zj's from the lower 
state is of course a somewhat arbitrary choice; it ensures the greatest accu- 
racy possible for the lower lying states (known as "state-specific" method), 
an alternative would be to choose the Zi so as to solve a least squares problem 
and minimize the overall error. 

To summarize, the minimal CI model for Li, Be, B, C, N, O, F, 
Ne excited states is as follows. Below, V^'^'^{Z_) denotes the symmetry 
subspace G l^(^) | L^^ = L{L + 1)^, S^^ = S{S + 1)^, m = p^}, 
where L is a non-negative integer, S a non-negative half-integer, and p = ±1. 

(A') (Choice of a parametrized, asymptotically exact family of subspaces) 
As in (A) 

(B') (Subspace eigenvalue problem) For each symmetry subspace V^'^'''(Z_) C V{Z_) 
Ef'(Z_) := eigenvalues of P{Z)HP{Z) on V^'^'P{Z), 

'■= corresponding orthonormal eigenstates, 
where P{Z) = orthogonal projector of L^((M^ x 12)^) onto V{Z) 

(C) (Variational parameter determination) For each symmetry subspace V^'^'P{Z^ C V{Z^ 
Z, := aigminzimm, Ef'{Z)), Ef := Ef'{Z,), := ^'^^(ZJ 

Let us summarize the additional properties of the model (A'), (B'), (C) 
beyond those of general symmetry-preserving CI (Lemmas 12.11 12.21) in a 
theorem. 

Theorem 2.1. Let N e {3, .., 10}, Z > N. The minimal CI model (A'), 
(B'), (C) has the following properties. 

(i) (Virial theorem) Any lowest normalized eigenstate of the model in 
a symmetry subspace (i.e., a joint eigenspace of the symmetry operators L^, 
S^, R) satisfies 

where T , V are the kinetic respectively potential part of the Hamiltonian ([^. 
(a) (Correct asymptotic behaviour) For N fixed and Z — oo, 

lim^ = l, lim ^ = 1, lim 1 1 iPf^ - R| 1 1 = 0, 

Ej AEj ' 

where E^^ < ■ ■ ■ < E^^ and Ei < ■ ■ ■ < E^ are the CI eigenvalues respec- 
tively the lowest eigenvalues of the Schrodinger equation (J\), AE'-''^ and AEj 
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are the spectral gaps E'^^ — E'^^ and Ej — Ei (j > 2), Pf^, Pj denote the 
projectors onto the corresponding eigenspaces, and \ \\-\\ \ is the operator norm 
on the N -electron Hilbert space Ll{{M.^ x ^2)^). 

Proof (i) follows from the fact that the manifold ^ze{o,ooy^y^'^'''^{Z_) is invari- 
ant under dilations \E' 1— » ^ x{xi, si, .., xn, sn) = A^^/^\l/(Axi, si, .., Axat, sat), 
A > 0, which makes the usual proof of the virial theorem applicable: normal- 
ized minimizers of (^aI-^^I^a) in this manifold satisfy = ^ (\E'a|-?^|^I/a)- 

(ii) is a consequence of the asymptotic results in |FG09j together with 
the elementary inequalities Ej < E*^^ < E^'^, where the E^'^ are the lowest 
eigenvalues of the PT model |FG09j . 

We remark that statement (ii) fails when the Slater orbitals f ill I) are replaced 
by finite linear combinations of Gaussians, or indeed by any functional form 
which fails to reproduce ( |TT|) asymptotically |FG09j . 

It is instructive to compare the above argument in favour of Slater or- 
bitals to the well known Kato cusp condition argument. Theorem 12.11 (ii) 
concerns the limit N and x general, Z ^ 00, whereas the asymptotic regime 
of the Kato cusp condition is N and Z general, x — *• 0; the latter is therefore 
insufficient to specify whole orbitals, as it only concerns their behaviour at 
x = 0. 

Finally, let us formulate a hierarchy of higher and higher dimensional 
CI models for the atom/ion with electrons which satisfy require- 
ments 1., 2., 3. The models are parametrized by the number rimax of included 
single-electron "shells", and the only modification compared to (A'), (B'), 
(C) is an enlargement of the family of subspaces V{Z_) in Step (A'), as fol- 
lows. For n = 1, ...,nmax, i = 0, ..,n-l, let Vne{Z_) := Span {i?„£(2'f, .., Z^,r) 
Xem{0, (j))6a{s) I m = —£, ..,£, a = ±|}. Here Z is the vector of dilation pa- 
rameters Z^, (r, (j), 6) are spherical polar coordinates in M^, the Rni{Zf, .., Z^, •) 
are orthonormal functions in L^(0, 00) with respect to the measure r'^dr which 
reduce to the usual radial hydrogen eigenfunctions when Z^ = ... = Z^ = Z , 
and the xem are spherical harmonics (see [FGOQj ). Then take 

ViZ) :=Span{|JC^}, 

d 

where d = {dn£)n=i,..,nmax,t=o,..,n runs over all partitions of A^, i.e. dni > 0, 
'^ni^rii = N. The minimal model (A'), (B'), (C) corresponds to taking 
nmax = 2 (i.e., including only the first and second "shell"), and imposing the 
additional condition that the number di^ of electrons in the subspace Vi^o 
equals two (i.e., assuming that the first shell is completely "filled"). 
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3 Minimal CI atomic energy levels and eigen- 
states 



3.1 Exact solution for given dilation parameters 

The key point allowing to solve the model (A'), (B'), (C) is the observation 
that the CI matrix P{Z_)HP{Z_) in a simultaneous eigenbasis of V{Z_) of 
the symmetry operators can be explicitly determined, and decouples into 
small invariant blocks. More precisely, as noted in |FG09] , exact expressions 
can be derived for the joint eigenstates of ([7j) and their matrix elements 
(\l/a|if|\l//3) in terms of one-body. Coulomb and exchange integrals of the one- 
electron orbitals (fTTll : and when restricting without loss of generality to 5*3 
maximal and L3 = the largest non-diagonal block is 2x2. For convenience 
we include the eigenfunctions and symbolic matrix elements in Tables [2H3] 
below. The symmetry type of the wavefunctions is also shown in Chemist's 
notation, which encodes the eigenvalues L{L + 1), S{S + 1) and p of L^, 
and R by the symbol '^^+'^X'^, where L corresponds to X via 0^5,. 1 ^ P, 
2 D, and no superscript u means p = 1, while u = o (for odd) stands for 
p = — 1. Recall the standard notation for one- and two-body integrals 

(a|6) = {a\hb), {ab\cd) = f dxidx2a* {xi)b{xi)- 7C*{x2)d{x2), (15) 

\Xi — X2\ 

where h is the one-body Hamiltonian — |A — Z/\x\. 

It remains to evaluate the one-body Coulomb and exchange integrals for 
the basis ( ITTl) . Despite the basis not being Gaussian, they can be evalu- 
ated exactly, by the method introduced in [FG09j : by Fourier calculus, we 
can re-write {ab\cd) = (27r^)~^/jg3 \k\~'^{ab*)*c*d; we then derive the Fourier 
transform of the pointwise products of the orbitals ffTTj) (see Lemma EU]), re- 
duce to ID integrals with the help of spherical polar coordinates in fc-space, 
and evaluate the remaining ID integrals - whose integrands turn out to be 
rational functions - via the residue theorem (or MAPLE). The result is as 
follows. 

Lemma 3.1. The Fourier transforms of pointwise products of the one-electron 
orbitals ( I77]j are as follows. In all cases j, £ = 1,2, 3, j 7^ £. 



Function 


Fourier Transform 


4>2s4>2a 
4>ls<l>2a 

't>2pj<t>2pj 
'i>la<t>2p, 
<i>2a<t>2pj 

4>2pj4>2pt 


wzf 

(4Zf + |fcp)2 

Z| f 2(Zi^2Z2) Z2(2Zi+Z2)(2Zi+5Z2) . 2Z|(2Zi+Z2)2 \ 


{4zi-2ZiZ2+z'^) \ {z'^+m^ (z'i+m^ ' {z'i+m" ) 

^/TiZl^^zy^ I 4(2Zi+Z2) (2Zi+Z2)3 \ 
^{4Zf-2ZiZ2+Z|) V3((Zi+Z2/2)2 + |fcP)2 3((Zi+Z2/2)2 + |*P)3 j 
Z| 6Z|fc| 

(Zj + |fcP)3 (Z| + |fcP)4 

2y2iZj'^zJ'^(2Zi+Z3)/cj 
((Zi+Z3/2)2 + |fc|2)3 
V3iZ-"'^zf''^ /8(Z2+Z3)2(2Zi+Z2)*i (32Zi+64Z2+48Z3)S:j \ 
16V4Z2-2ZiZ2 + Z| V (((Z2 + Z3)/2)2 + |)=P)'l 3(((Z2 + Z3)/2)2 + |fcP)-' j 

(Z| + |fc|2)l 
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R 








Li 





3 
4 


1 




|1T2> 


2(1|1) + (2|2) + (11|11) +2(11|22) - (12|21) 




2 


1 


-1 




|1T3> 


2(1|1) + (3|3) + (11|11) +2(11|33) - (13|31) 


Be 





n 


1 


'S 


|1T22) 


2(1|1) + 2(2|2) + (11|11) + 4{11|22) - 2(12|21) + (22|22) 












^(|1T33) + |1T44) + |1T55)) 


2(1|1) + 2(3|3) + (11|11) + 4(11|33) - 2(13|31) + (33|33) + 2(34|43) 














\/3(23|32) 




2 


n 


_1 


1 po 


— ('llT23\ — nT23\l 


2(1 1) + (2|2) + (3 3) + (11 11) + 2(11 22) — (12 21) + 2(11 33) 
-(13|31) + (22|33) + (23|32) 




2 


2 


-1 


3po 


|1T23) 


2(1|1) + (2|2) + (3|3) + (11|11) + 2(11|22) - (12|21) + 2(11|33) 
-(13|31) + (22|33) - (23|32) 








1 


3p 


|1T45> 


2(1|1) + 2(3|3) + (11|11) + 4(11|33) - 2(13|31) + (33|44) - (34|43) 




6 





1 




i(2|lT33>-|lT44)-|lT55)) 


2(1|1) + 2(3|3) + (11|11) + 4(11|33) - 2(13|31) + (33|33) - (34|43) 


B 





4 


1 




^ (|lT233) + 1T244) + |lT255)) 


2(1 1) + (2|2) + 2(3|3) + (11|11) + 2(11|22) - (12|21) + 4(11|33) 
-2(13|31) + 2(22|33) - (23|32) + (33|33) + 2(34|43) 




n 


T 


-1 


Ago 


iTsir,) 


2(1|1) + 3(3|3) + (11|11) + 6{11|33) - 3(13|31) + 3(33|44) - 3(34|43) 




2 


3 
4 


-1 


2po 


1T223) 


2(1|1) + 2(2|2) + (3|3) + (11|11) +4(11|22) - 2(12|21) + 2(11|33) 
-(13|31) + (22j22) + 2(22|33) - (23|32) 












^(|1T344> + |1T355)) 


2(1|1) + 3(313) + (lljll) + 6(11|33) - 3(13|31) + (33|33) + 2(33|44) 














\/2(23|32) 








1 


^-P 


^ (2|11245) - |lT245) - |lT245» 


2(1|1) + (2|2) + 2(3|3) + (11|11) + 2(11|22) - (12121) + 4(11|33) 
— 2(13|31) + 2(22|33) + (23 32) + (33|44) — (34|43) 




2 


15 


1 


JP 


1T245) 


2(1|1) + (2|2) + 2(3|3) + (11|11) + 2(11|22) - (12|21) + 4(11|33) 
-2(13|31) + 2(22|33) - 2(23|32) + (33|44) - (34|43) 




6 


3 


1 




^ (2|11233) - |11244) - |11255» 


2(1|1) + (2|2) + 2(3|3) + (11|11) + 2(11|22) - (12|21) + 4(11|33) 
-2(13|31) + 2(22|33) - (23|32) + (33|33) - (34|43) 








-1 




^ (2|11345) - |1T345> - |lT345)) 


2(1|1) + 3(3|3) + (11|11) + 6(11|33) - 3(13|31) + 3(33|44) 


C 








1 




^ (|lT2233> + |lT2244) + |lT2255» 


2(1|1) + 2(2|2) + 2(3|3)0 + (11|11) + 4(11|22) - 2(12|21) + 4(11|33) 
-2(13|31) + (22|22) + 4(22|33) - 2(23|32) + (33|33) + 2(34|43) 












^ (|1T3344) + |1T3355) + |lT4455» 


2(1|1) + 4(3|3) + (11|11) + 8(11|33) - 4(13|31) + 2(33|33) + 4(33|44) 














2(23|32) 




n 


2 


-1 




^ (3 1 112345) - |1T2345) 


2(1|1) + (2|2) + 3(3|3) + (11|11) + 2(11|22) - (12|21) + 6(11|33) 
— 3n 3131 ^ 4- 3f99l33i + f23l321 A- 3f 331441 — 3^341431 







6 


-1 




|1T2345> 


2(1|1} + (2|2) + 3(3|3) + (11|11) + 2{11|22) - (12|21) + 6(11|33) 




2 





-1 


ipo 


|(|1T2344> - 1112344) 

+ lizooo/ — \lLZooO)) 


2(1|1) + (2|2) + 3(3|3) + (11|11) + 2(11|22) - (12|21) + 6(11|33) 
— o(io|oij + o\2.L\oo) -r \6o\oo) -r zloo 44J 




2 


2 


1 


1 


1 1 "Too" /I rz \ 


2(l l) + 2(2 z)U + •^V'^PJ "T (l-l-|llj + 4(li|z2j — 2(lz|zl) +4(ll|3o) 
-2(13|31) + (22|22) + 4(22|33) - 2{23|32) + (33|44) - {34|43) 














9l'lll^-l-4('3l3l + f1lllll4-Sflll33l — 4ri3l31'\-l- ('331331 4- 5('33l44l 
-3(34|43) 












cross 


(23|32) 








-1 


3po 


^(|1T2345) + |1T2355)) 


2(1|1) + (2|2) + 3(3|3) + (11|11) + 2(11|22) - (12121) + 6(11|33) 
-3(13|31) + 3(22|33) - 2(23|32) + (33|33) + 2(33|44) 




6 





1 


'D 


^ (2|lT2233) - |lT2244) - 1112255)) 


2(ljl) + 2(2|2) +2(3|3) + (lljll) + 4(11|22) - 2(12|21) + 4(ll|v33) 
-2(13|31) + (22|22) + 4(22|33) - 2(23|32) + (33|33) - (34143) 












^ (21114455) - I1T3344) - IIT3355)) 


2(1|1) + 4(3|3) + (11|11) + 8(11|33) - 4(13|31) + 2(33|33) + 4(33|44) 
-3(34|43) 












cross 


-(23|32) 








-1 




^(2|lT2345) - IIT2345) - IIT2345) 
+2|112345) - |lT2345> - |lT2345)) 


2(1|1) + (2|2) + 3(3|3) + (11|11) + 2(11|22) - (12|21) + 6(11|33) 
-3(13|31) + 3(22|33) + 3(33|44) 




6 


2 


-1 




Jj(2|112345> - 1112345) - |112345)) 


2(ljl) + (2|2) + 3(3|3) + (11|11) + 2(11|22) - (12|21) + 6(11|33) 
-3(13|31) + 3(22133) - 2(23132) + 3(33144) 



Table 2: H matrix element expressions for Li-C in subspace with 5*3 maximal, 
L3 = 0; 'cross' denotes the off-diagonal term in the 2x2 matrix. Orbital 
notation is 1 := 01^, 2 := 02s, 3 := 02p3, 4 := 02pi, 5 := 02p2- Integrals (a|a), 
{ah\cd) are as in (fT^ . 
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L2 


s- 








(H) 


N 





4 


1 




^(|1T23344) + |1T23355> 
+|lT24455)) 


2(1|1) + (2|2) + 4(3|3) + (11|11) + 2(11|22) - (12|21) + 8(11|33) 
-4(13|31) + 4(22|33) - 2(23|32) + 2(33|33) + 4(33|44) 







15 
4 


-1 


4go 


|lT22345) 


2(1|1) + 2(2|2) + 3(3|3) + (11|11) + 4(11|22) - 2(12|21) + 6(11|33) 
-3(13|31) + (22|22) + 6(22|33) - 3(23|32) + 3(33|44) - 3(34|43) 




2 


3 
4 


-1 


2po 


^ (I1T22344) + I1T22355)) 


2(1|1) + 2(2|2) + 3(3|3) + (11|11) + 4(11|22) - 2(12|21) + 6(11|33) 
-3(13131) + (22122) + 6(22133) - 3(23132) + (33133) + 2(33144) 












|lT34455> 


2(1|1) + 5(3|3) + (11|11) + 10(11|33) - 5(13|31) + 2(33|33) 
+8(33|44) - 4(34|43) 












cross 


\/2(23|32) 








1 


'^P 


^(2| 1123345) - |lT23345) 
-|1T23345)) 


2(1|1) + (2|2) + 4(3|3) + (11|11) + 2(11|22) - (12|21) + 8(11|33) 
-4(13|31) + 4(22|33) + (33|33) + 5(33|44) - 3(34|43) 




2 


15 
4 


1 


4 P 


|liz334o; 


2(l|Ij + (2|z) + 4^313) + (li|lij + 2(li|z2J — (^iz|2ij + 8(ii|33; 

-4(13|31) + 4(22|33) - 3(23|32) + (33|33) + 5(33 44) - 3(34 43) 




fi 


4 


-1 


■Ijjo 


^(2|lT22345> - |lT22345) 
-|1T22345)) 


2(1|1) + 2(2|2) + 3(3|3) + (11|11) + 4(11|22) - 2(12|21) + 6(11|33) 
-3(13|31) + (22|22) + 6(22|33) - 3(23|32) + 3(33|44) 








1 




-i5(2|lT24455) - |lT23344) 
-|1T23355)) 


2(1|1) + (2|2) + 4(3|3) + (11|11) + 2(11|22) - (12|21) + 8(11|33) 
-4(13|31) + 4(22|33) - 2(23|32) + 2(33|33) + 4(33|44) - 3(34|43) 











1 


'S 


^(1 11223344) + |1T223355) 
+I1T224455)) 


2(1|1) + 2(2|2) + 4(3|3) + (11|11) + 4(11|22) - 2(12|21) + 8(11|33) 
-4(13|31) + (22|22) + 8(22|33) - 4(23|32) + 2(33|33) + 4(33|44) 












|1T334455) 


2(1|1) + 5(3|3) + (11|11) + 12(11|33) - 6(13|31) + 3(33|33)) 
+ 12(33|44- 6(34|43) 












cross 


^(23|32) 




2 





-1 


ipa 


(|lT234455) - |lT234455)) 


2(1|1) + (2|2) + 5(3|3) + (11|11) + 2(11|22) - (12|21) + 10(11|33) 
-5(13|31) + 5(22|33) - {23|32) + 2(33|33) + 8(33|44) - 4(34|43) 




2 


2 


-1 




|lT234455) 


2(1|1) + (2|2) + 5(3|3) + (11|11) + 2(11|22) - (12|21) + 10(11|33) 
-5(13|31) + 5(22|33) - 3(23|32) + 2(33|33) + 8(33|44) - 4(34|43) 










3p 


|iT99'->^4c;\ 


-4(13|31) + {22|22) + 8{22|33) - 4{23|32) + (33|33) + 5(33|44) 
-3(34 43) 




(i 


{) 


]^ 


1 p) 


^^2nT9944fi^\ — hT923344\ 
-|1T223355)) 


-4(13|31) + (22|22) + 8(22|33) - 4(23|32) + 2(33|33) + 4(33|44) 
-3(34 43) 


F 





3 
4 


1 




|lT2334455) 


2(1|1) + (2|2) + 6(3|3) + (11|11) + 2(11|22) - (12|21) + 12(11|33) 
-6(13|31) + 6(22|33) - 3(23|32) + 3(33|33) + 12(33|44) - 6(34|43) 




2 


3 
4 


-1 




1 1T2234455) 


2(1|1) + 2(2|2) + 5(3|3) + (11|11) + 4(11|22) - 2(12|21) + 10(11|33) 
-5(13|31) + (22|22) + 10(22|33) - 5(23|32) + 2(33|33) + 8(33|44) 
-4(34 43) 


Nc 








1 


IS 


|1T22334455) 


2(1|1) + 2(2|2) + 6(3|3) + (11|11) + 4(11|22) - 2(12|21) + 12(11|33) 
-6(13|31) + (22|22) + 12(22|33) - 6(23|32) + 3{33|33) + 12(33|44) 
-6(34 43) 



Table 3: H matrix element expressions for N-Ne in subspace with 5*3 maxi- 
mal, L3 = 0; 'cross' denotes the off-diagonal term in the 2x2 matrix. Orbital 
notation is 1 := 01^, 2 := 02s, 3 := (^2^3, 4 := (^2pi, 5 := (^2^2- Integrals (a|a), 
{ab\cd) are as in (JT5l) . 
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Lemma 3.2. Using the abbreviated notation 1 = (pis, 2 = 02s; 3 = 4'2pz, 
4 = 02pi; the one-body, Coulomb and exchange integrals ( 173]) for the orbitals 
07]) respectively the PT orbitals (Zi = Z2 = Z'^ = Z ) are given by 



Integral 



m 

(2|2) 
(3|3) 



(11|11) 
(11|22) 
(12|21) 
(22|22) 

(11|33) 
(13|31) 
(22|33) 

(23|32) 

(33|33) 
(33|44) 
(34|43) 



CI value 



\Z\ — ZZi 



4Zf-2ZiZ2+7Zj _ ZZ-2 iZ{-iZiZ2+3Z^ 
24 42^-221^2+21 4 iZ'f-2ZiZ2+Z'^ 

1^1 ~ i^^3 



2i 22(82? +42? 22+42i2|+Z|) 
{22i +22)3 (42J -22i 22 +2|) 

162i'2| 

(42J -22i 22 +2|) (22i +22 
22 (14882;'-19522j' 22+17522? 2f-8402i2j+2452|) 
512 {42?-22i22+2|)2 

2i23(82?+202? 23+202? 2f+102i2|+2|) 
(22i+23)= 
1122? 2| 
3(22i+23)' 

(4(42?-22i2^+%)(22+23)^) K"^? " ^^'^^ + ^^'X^^' + '^^^^a + 21Z|2| + 352|2|) 
+32|2|(282? - 282i22 + 1122) ^jZ2Zl(4Z^ - 42i22 + 32|) + 2|{42? - 42i22 + 32|))] 
2| 2| (7402? + 1522i 22 + 172| -4222 23 - 5882i 23 + 1262| ) 
(9(22 +23)''(42?-22i22+2|)) 
^Z, 



PT value 



16 7 
729 

77 7 

512 

59 y 
243-^ 
112 7 
6561 



83 7 
512-^ 
15 7 

501 7 
2560 
447 7 
2560 
27 7 
2560 



This table together with Tables [2H3] yields, for any given values of the Z,, 
the exact solution of the linear part (B') of the CI model in the nondegenerate 
symmetry subspaces. 

In the 2D subspaces, the tables need to be combined with the analytic 
expression for the eigenvalues of the 2x2 matrices {{^ i\H\^ j))^^^^^ (see 
jphnO] and denote {^i\H\^<j) =: Hij), 



A. 



22 



± 



(16) 



and corresponding normalized eigenstates. 



■ (17) 



12 



Thus we have analytic expressions for all eigenvalues and eigenvectors of 
P{Z)HP{Z) in terms of the Zi. 



3.2 Numerical optimization of dilation parameters 

The final stage is to minimize the exact energy levels over the Zj (Step (C) 
of the minimal CI model), which is performed using MAPLE. Since we are 
dealing with only a 3-parameter minimization over explicit rational or square 
root functions, we obtain highly accurate numerical energy levels, along with 
their eigenspaces and symmetries. In particular, all digits indicated in the 
Tables below are beheved to be exact relative to the underlying model (A'), 
(B'), (C). 



17 



3.3 Final result 



The minimal CI energy levels, along with the minimizing values of the di- 
lation parameters Zj, for N = 3, . . . , 10, are shown in Tables H] and The 
corresponding eigenspaces are as given in Tables [2HS1 

4 Comparison with large-scale numerical cal- 
culations and experiment 

4.1 Ground state energies and ground states 

The results in Tables H] and [5] show that the symmetry of the ground state of 
the model (A), (B), (C) agrees with experiment in every case, and that the 
ground state energies capture around 99 % of the experimental energy. 

We consider this agreement very good for such a low-dimensional projec- 
tion of the Schrodinger equation. In the case of Beryllium, our ground state 
CI energy even outperforms the benchmark numerical multi-determinant 
Hartree-Fock results of |TTST94] . This demonstrates that a careful choice 
of basis and considering the full Hamiltonian, including all correlation terms, 
can be more effective than large numerical computations. 

It is also of theoretical interest to compare with the best numerical values 
in the literature, which rely on more high-powered approaches. The table 
below compares, in a typical example, our asymptotics-based minimal CI 
results, the MDHF results of Tatewaki et al. (also based on a small number 
of determinants but on a huge one-electron basis set, considered essentially 
complete), the MPII results of Canal Neto, Muniz, Centoducatte and Jorge, 
and the benchmark Full CI results of Bauschlicher and Taylor. 

Other examples we considered gave a similar picture. In particular, for 

> 5 asymptotics-based minimal CI was never outperformed by more than 
one digit in all tested cases, not even by the recent explicitly correlated, multi- 
configurational variational Monte Carlo results [GBS02J; for N = 4 (Be) the 
sophisticated iterative subspace recursions of ^BM04^ iNNKIOTj - which lead 
to complicated final wavefunctions with 3 x 10^ respectively 2 x 10^ DOF's 
- only yield energies which are one respectively two digits more accurate. 
While from an applications point of view an accuracy gain of one digit can 
be very important, the fact remains that the required computational effort 
is larger by many orders of magnitude. A tentative conclusion is that a 
significant part of the quality of quantum chemistry models lies in mak- 
ing a sophisticated initial ansatz, while subsequent efforts to include more 
and more contributions appear to exhibit the same disappointing scaling be- 
haviour expected from a direct discretization of a problem suffering from the 
curse of dimension. 

Also of theoretical interest is the large gain in accuracy of minimal CI over 
the PT model (i.e. first order perturbation theory with respect to electron 
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Z3 
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Eex'p 


Emdhf 


EpT 


AEci 


^Eexp 


^Emdhf 


N 




-54.1597 


6.7117 


4.7535 


3.7924 




-54.6117 


-54.4009 


-49.1503 












-54.0407 


6.7124 


4.7711 


3.7317 




-54.5241 


-54.2962 


-48.9288 


0.1190 


0.0876 


0.1048 




2po 


-54.0075 


6.7110 


4.7893 


3.7162 


-0.2091 


-54.4803 


-54.2281 


-48.8195 


0.1523 


0.1314 


0.1728 




ip 


-53.7666 


6.6854 


4.8658 


3.7592 




(-54.2101) 


-53.9883 


-48.1630 


0.3932 


(0.4016) 
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4.8414 


3.7065 
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3.6669 
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4.7893 


3.7162 
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5.5613 
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-74.8094 


-66.7048 












-74.3004 


7.7122 


5.5709 


4.3828 




-75.0357 


-74.7293 


-66.5360 


0.0928 


0.0723 


0.0801 






-74.2328 


7.7103 


5.5967 


4.3628 


-0.2283 


-74.9540 


-74.6110 


-66.3421 


0.1603 


0.1540 


0.1984 




^ po 


-73.7784 


7.6805 


5.6490 


4.3916 




(-74.5324) 


-74.1839 


-65.3265 


0.6147 


(0.5756) 


0.6255 
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-73.4204 


7.6785 


5.5620 


4.3549 






-73.8720 


-64.8578 


0.9727 




0.9374 






-72.8054 


7.7103 


5.5967 


4.3628 


4.3811 






-63.4984 


1.5877 
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6.3576 


5.0587 




-99.8060 


-99.4093 


-87.6660 
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-97.8704 


8.6748 


6.4189 


5.0416 




(-99.0322) 


-98.5312 


-85.8342 


0.8800 


(0.7738) 


0.8781 
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-127.5695 


9.7101 


7.1469 


5.7177 
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Method 


1st order PT \\<VA)i')\ 


iniiiimal CI (this paper) 


MDHl-' |'l"i'H'r!)4| 


MPll K!NM(!.l(iK| 


FCI PTO)| 


DOF's 


8 


11 


~ 10^ {estimate) 


not given 


2.8x10" 


Error 


12% 


1.06% 


0.40% 


0.28% 


0.21% 



Table 6: Fluorine atom ground state energy (various methods). DOF's means 
the number of computational degrees of freedom. 



interaction) [FG09], since the two models differ only by the optimization step 
(C) over dilation parameters. 



Atom 


Li 


Be 


B 


C 


N 





F 


Ne 


PT Error 
CI Error 


5.6% 
0.9% 


6.2% 
0.6% 


7.8% 
0.7% 


9.0% 
0.8% 


10.0% 
0.8% 


11.2% 
1.0% 


12.2% 
1.0% 


13.0% 
1.1% 



Table 7: Percentage error in the ground state energy for the PT and minimal 
asymptotics-based CI method. 

Some insight can be gained from comparing the CI orbitals resulting from 
energy minimization with the "bare" PT orbitals. It is clear from Tables 
m and that Zi ^ Z — 0.3 and hence the PT model Is orbitals are a fair 
approximation to those in the CI model. But this is not true for the 2s and 
2p orbitals since Z2 is lower than Z by about 2, and Z3 is lower by about 2 
to 4. 

Physically this is intuitive from the idea that the Is orbitals partially 
screen the nuclear charge felt by the 2s and 2p orbitals, making the 2s and 
2p electrons behave as they would in the potential of a nucleus with reduced 
nuclear charge. 

Mathematically, one can at least explain why the Zi differ from their PT 
value of Z. The CI wavefunctions satisfy the virial theorem (see Section [2]); 
by contrast the deviation of the PT wavefunctions from the correct virial 
ratio between potential to kinetic energy of —2 is large, because these states, 
being ground states of a non-interacting Hamiltonian, have a ratio of —2 
for (potential energy without electron repulsion) to kinetic energy. (From 
|FG09j . the actual virial ratios of the PT ground states for Li, Be, B, C, 
N, O, F, Ne are -1.6969, -1.6881, -1.6615, -1.6379, -1.6173, -1.5956, -1.5778, 
-1.5615.) 

We now discuss the obtained wavefunctions. Our work provides for the first 
time few-parameter, explicit, closed-form wavefunctions for the low- lying 
eigenstates of the atoms Li, Be, B, C, N, O, F, Ne which are of chemically 
relevant accuracy. These can be used as a source of numerous theoretical 
insights. 

As an important application, the wavefunctions given by Tables |2H3] and eq. 
( |T71) and their ordering given in Tables HHH] confirm and make quantitative 
the qualitative mathematical picture of the periodic table obtained in [FG09j 
by asymptotic analysis of the Schrodinger equation for strongly positive ions. 
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For instance, they affirm the conclusion of |FG09] that the empirical shell 
ordering rule of quantum chemistry (as the electron number increases, the 2s 
shell is "ffiled" before the 2p shell) is only correct in a probabilistic sense. In 
degenerate symmetry subspaces, the minimal CI eigenstates contain the two 
configurations ls^2s^2p^~'^ and ls^2p^~^ (see the discussion of 2s-2p reso- 
nance in |FG09] ). The state with lower energy is dominated by the ffist con- 
figuration, i.e. the coefficient for the part of the wavefunction in ls^2s^2p^~^ 
is larger than the part in ls^2p^~^. The reverse is true for the higher energy 
state. Nevertheless, the minority contributions are of significant size (36%, 
27% and 17% in case of the B, Be, C ground state). 

4.2 Spectral gaps and ionization energies 

These are an extremely tough test of any model, due to the multiscale effect 
that they are smaller by two to three orders of magnitude (see Table [T]). 

First, note how our eigenstate tables allow to trace spectral gaps to the 
size of individual Coulomb and exchange integrals, revealing the cancellations 
that lead to the small size of gaps compared to total energies (see Table [1]). 

As an example of a 2s-2p spectral gap, consider the ^5* ground state and 

first excited state of Lithium. Table [2] shows that the gap at fixed values 
of Zi, ^2, is given by the difference in one-body energy and interaction 
with the Is shell of the 2p and 2s orbitals, [(3|3)+2(11|33)-(13|31)]-[(2|2) + 
2(11|22) — (12|21)]. Substituting for simplicity the bare values Zi = Z = "i 
into the table in Lemma [3^ the difference between the Coulomb terms is only 

= 0.19753... (and that between the exchange terms only = 0.01463...), 
which is much smaller than the common part 2(1|1) + (2|2) + (11|11) = 
2(11) + (3|3) + (11|11) = -f = -8.25 contained in each of the states. 

As an example of an energy level splitting between two states with an 
equal number of Is, 2s and 2p orbitals, consider the "^S" ground state and 
ffist excited state of Nitrogen. A look at Table [3] reveals that the energy 
difference consists only of the exchange term — 3(34|43), which is present in 
the ground state due to the parallel spins of the three p-orbitals, but absent 
in the excited state. 

Next, as shown in Tables H] and [5], the spectral gaps for the CI model 
are in good agreement with experimental data (most are within 10%) and 
comparable to the predictions of numerical studies with a much larger number 
of degrees of freedom |TTST94j . Considering for example the first three 
spectral gaps of Nitrogen, Carbon and Oxygen, CI has the more accurate 
value in five out of nine cases, and the less accurate value in the remaining 
four cases. 

To achieve this accuracy, the minimal form (C) of relaxation of orbitals in 
the CI model is needed, as the "bare" PT orbitals, despite sharing asymptotic 
exactness in the large nuclear charge limit, give very poor spectral gaps, with 
errors in the order of 200 — 300%. 

One interesting qualitatively new feature introduced by the CI model 
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10 15 20 . 25 



Figure 2: Minimal CI energy differences -£"350 — Eijjo for the Carbon isoelec- 
tronic sequence N = 6, Z = 6,. ..,28. The predicted existence of a level 
crossing is confirmed by experiment |RJK"'"07] . 

is the possibility for energy levels to cross as the nuclear charge Z varies 
(see Figure [2]) . This is due to the non-linearity of the energy levels in Z 
arising from the minimization over the dilation parameters Z_ = {Zi, Z2, Z^). 
(Note that for Zi = Z2 = Z^ = Z, the energy levels have the special form 
—aZ"^ + biZ |FG09] . yielding linearity of gaps in Z.) This enables us to 
discuss, for example, the and states of the Carbon isoelectronic 
sequence. We recall from |FG09j that both Hund's rules and the Hartree-Fock 
picture predict the universal ordering Esgo < Eijjo, which agrees with the 
experimental orderings for Carbon. However, for Z > 20 the experimental 
ordering is found to be reversed. This crossing is beautifully captured by the 
minimal CI model, this time for Z > 23. 

We now discuss another important class of energy differences, ionization 
energies. The latter are defined to be (writing Ei{N, Z) to indicate the de- 
pendence of the ground state energy on the number of electrons and the 
nuclear charge) I{N, Z) = Ei{N — 1, Z) — Ei{N, Z). Physically this corre- 
sponds to the energy required to remove one electron from a system with 
nuclear charge Z and N electrons. The calculated first ionization energies 
I{N, N) of the minimal CI Model, in atomic units, are as follows: He 0.8477, 
Li 0.1912, Be 0.3237, B 0.2346, C 0.3142, N 0.3960, O 0.2708, F 0.3958, 
Ne 0.4141. The experimental ionization energies i Huh93] are: He 0.9036, Li 
0.1980, Be 0.3426, B 0.3049, C 0.4138, N 0.5341, O 0.5000, F 0.6402, Ne 
0.7925. 

Figure [S] shows that the qualitative prediction for the ionization energies 
is very good when compared to experimental data. In particular, all local 
minimizers (H, Li, B, O), local maximizers (He, Be, N), global minimizers 
(Li) and global maximizers (He) are predicted correctly. This is all the more 
remarkable when remembering that tiny eigenvalue differences for partial 
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H He Li Be B C N O F 

Atom 



Figure 3: Experimental (dashed) and minimal CI (solid) ionization energies. 

differential operators on very high-dimensional spaces up to M^^ are under 
consideration here. 

Quantitatively, for the smaller atoms our results are comparable to (and 
in case of Be better than) MDHF calculations with much larger basis sets up 
to 20sllp |JAH01j . For the larger atoms the minimal dimensionality of our 
CI subspace finally makes itself felt, and a larger subspace (e.g. as described 
at the end of Section 2.4) would be needed to make the qualitative agreement 
quantitative. 

Again, it is also instructive to compare with the PT model [FG09j . Its 



ionization energies, which are easily read off from the exact results of |FG09j , 
even turn out to have the wrong sign. This shows that relaxation of orbitals 
is important for the description of ionization processes, and that the relax- 
ation step (C) in the minimal asymptotics-based CI model is essential for 
understanding the nontrivial experimental graph in Figure [3l 
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